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We study the quantum effects induced by bulk scalar fields in a model with a de Sitter (dS) brane 
in a flat bulk (the Vilenkin-Ipser-Sikivie model) in more than four dimensions. In ordinary dS space, 
it is well known that the stress tensor in the dS invariant vacuum for an effectively massless scalar 
(?Tie ff = m 2 + £1Z — with 1Z the Ricci scalar) is infrared divergent except for the minimally coupled 
case. The usual procedure to tame this divergence is to replace the dS invariant vacuum by the Allen 
Follaci (AF) vacuum. The resulting stress tensor breaks dS symmetry but is regular. Similarly, in 
the brane world context, we find that the dS invariant vacuum generates (T MI/ ) divergent everywhere 
when the lowest lying mode becomes massless except for massless minimal coupling case. A simple 
extension of the AF vacuum to the present case avoids this global divergence, but (T M „) remains 
to be divergent along a timelike axis in the bulk. In this case, singularities also appear along the 
light cone emanating from the origin in the bulk, although they are so mild that {T^ v ) stays finite 
except for non-minimal coupling cases in four or six dimensions. We discuss implications of these 
results for bulk inflaton models. We also study the evolution of the field perturbations in dS brane 
world. We find that perturbations grow linearly with time on the brane, as in the case of ordinary 
dS space. In the bulk, they are asymptotically bounded. 

YITP-04-38, KUNS-1927 

1. INTRODUCTION 

The brane world (BW) scenario 0,0 has been intensively studied in the recent years. Little is known yet concerning 
the quantum effects from bulk fields in cosmological models 0- [3- S El ■ Quite generically, one expects that local 
quantities like (T^„) or ((f> 2 ) can be large close to the branes, due to the well known divergences appearing in Casimir 
energy density computations. This has been confirmed for example in for flat branes. These divergences are of 

ultraviolet (UV) nature and do not contribute to the force. Hence, they are ignored in Casimir force computations. 
However, they are relevant to the BW scenario since they may induce large backreaction, and are worth of investigation. 

In this article, we shall shed light on another aspect of objects like (T^) in BW. We shall point out that they 
can suffer from infrared (IR) divergences as well. These divergences arise when there is a zero mode in the spectrum 
of bulk fields in brane models of RSII type with dS branejl], |||. The situation is analogous to the case in dS space 
without brane. It is well known that light scalars in dS develop an IR divergence in the dS invariant vacuum. The 
main purpose of this article is to explore the effects of scalar fields with light modes in a BW cosmological setup of 
the RSII type Q. To consider massless limit of scalar field in inflating BW is especially well motivated in the context 
of 'bulk inflaton' models 0, IS H> 0> 0] i in which the dynamics of a bulk scalar drives inflation on the brane. In the 
simplest realizations, the brane geometry is close to dS and the bulk scalar is nearly massless. 

Let us recall what happens in the usual dS case ^3|. For light scalars m e s <C H (with H the Hubble constant) in 
dS, (<f> 2 ) and {T^ v ) in the dS invariant vacuum develop a global IR divergence ~ l/m 2 ff . To be precise, this depends 
on whether the field is minimally coupled or not. What we have in mind is a generic situation in which the effective 
mass m 2 ff = m 2 +£1Z is small, and £ ^ 0. In these cases (T^ v ) diverges as mentioned. The point is that in the generic 
massless limit, another vacuum must be chosen to avoid the global IR divergence. This process breaks dS invariance 
[l4|. but this shall not really bother us. The simplest choice is the Allen Follaci (AF) vacuum , in which the stress 
tensor is globally finite and everywhere regular. The massless minimally coupled case is special [l5|| , and it accepts a 
different treatment which gives finite (T^ u ) without violating dS invariance. 

In the BW scenario pj, the bulk scalar is decomposed into a continuum of KK modes and bound states. Here we 
consider the case that there is a unique bound state with mass mj. If is light, (</> 2 ) and (T^) for the dS invariant 
vacuum will also diverge like \ jmh. In this case, again, one will be forced to take another vacuum state like the AF 
vacuum. Then one naive question is what is the behavior of the stress tensor in such a vacuum in the BW. Also, one 
might expect singularities on the light cone emanating from the center (the fixed point under the action of dS group) 
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if we recall that the field perturbations for a massless scalar in dS grow like (<f> 2 ) ~ \, where \ is the proper time in 
dS [SHI 13 E3 (see also |H). The light cone in the RSII model corresponds to x — * °°- 

Before we start our discussion, we should mention previous calculation given in Ref. [2]]]. In that paper the stress 
tensor for a massless minimally coupled scalar was obtained in four dimensions, in the context of open inflation. 
Montes showed that (T^) can be regular everywhere except on the bubble. As we will see, these properties hold as 
well in other dimensions, but only for massless minimal coupling fields. 

For simplicity, we consider one extremal case of the RSII model in which the bulk curvature and hence the bulk 
cosmological constant is negligible. We take into account the gravitational field of the brane by imposing Israel's 
matching conditions. The resulting spacetime can be constructed by the 'cut-and-paste' technique. Imposing mirror 
symmetry, one cuts the interior of a dS brane in Minkowski and pastes it to a copy of itself (see Fig.QJ. Such a model 
was introduced in the context of bubble nucleation by Vilenkin p2^ | and by Ipser and Sikivie |23| , and we shall refer 
to it as 'the VIS model'. 

This article is organized as follows. In Section we describe the VIS model and introduce a bulk scalar field with 
generic bulk and brane couplings. The Green's function is obtained first for the case when the bound state is massive, 
rrid > 0. The form of (T^) in the limit m,d — > is also obtained. In Sectional we consider an exactly massless bound 
state TOrf = 0, and we present the divergences of the AF vacuum. The case when the bulk mass vanishes is technically 
simpler and explicit expressions for (T^ v ) can be obtained. This is done in Section With this, we describe the 
evolution of the field perturbations in Section [5J and conclude in Section |SJ 




FIG. 1: (a) Conformal diagram of the Minkowski space. The surfaces corresponding to constant r, \, T an d ip coordinates are 
indicated. The light cone r — corresponds to the dashed lines, (b) Conformal diagram of the VIS model. 



2. SCALAR FIELDS IN THE VIS MODEL 



In this Section we consider a generic scalar field propagating in the VIS model, describing a gravitating brane in 
an otherwise flat space • Specifically, the space time consists of two copies of the interior of the brane glued at 

the brane location, as illustrated in Fig. ^ 

In the usual Minkowski spherical coordinates the metric is ds 2 = —dT 2 + dR 2 + R 2 dfl 2 n y where dVl 2 n ^ stands for 
the line element on a unit n sphere. From the symmetry, it is convenient to introduce another set of coordinates. The 
Rindlcr coordinates, defined by 



R = r cosh x , 
T = r sinhx , 
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cover the exterior of the light cone emanating from R = T = 0. In terms of them the brane location is simply r = Tq, 
and the metric looks like 



ds 2 = dr 2 + r 2 dS 2 



(n+l) 



(1) 



where dS 2 n+1 ^ is the line element of de Sitter (dS) space of unit curvature radius. Thus, the Hubble constant on the 
brane is H = 1/ro. In order to cover the interior of the light cone, we introduce the 'Milne' coordinates according to 

R = t sinh %j) , 
T = t cosh ip . 



In these coordinates, the metric is ds 2 = —dr 2 + T 2 [dip 2 + sinh 2 ip dQ 2 n ->]. Note that we can go from the Rindler to 
the expanding (contracting) Milne regions making the continuation r = ±ir and \ — V' T (±)i7r/2. Here upper and 
lower signatures correspond to +ie and —it prescriptions, respectively. 

We consider a scalar field even under Zi symmetry, with generic couplings described by the action 



S = -\l [(d^f + (M 2 + £R) 2 ] - / 2£tr K] 4> 2 

1 ./Bulk ./Brane 



(2) 



where M and [i are the bulk and brane masses, 7Z is the Ricci scalar and tr/C is the trace of the extrinsic curvature. 
The latter arises because the curvature scalar contains S function contributions on the brane, and the factor 2 in front 
of it is due to the symmetry. The stress tensor for a classical field configuration can also be split into bulk and a 
surface parts as[l3l |24j 



T, 



ibulk 
brane 



= {i-2i) d^d v <t> 



l-4£ 



(d(j>f + (M 2 + £R,) f gilv - 2£R^ 2 - 2^V M V, 



[(4£ - 1) Meff hij - 2£ Kij\ S(r - r ) , 



(3) 



where hij is the induced metric on the brane, TZ^ the Ricci tensor, and the equation of motion has been used. 1 Here 
we have introduced an effective brane mass as 



fi eS = fi + 2(n + 



(4) 



where H = 1/ro is the Hubble constant on the brane. Then, the v.e.v. {T^ u ) in point splitting regularization is 
computed as 2 



jrp \bulk 



It 

2 V 



G {l) (x,x') 



with 



and 



Urn |(1 - 2£)d^' v - 9tlv {g^d x d' a + M 2 ) - 2£V M V„j , 



. bra 



(4C-l) Meff ^ i -2^y 5(r-r )G^(x,x')\. , 



(5) 



(6) 



(7) 



where d'^ = d/dx'^. This expression is extended to the case with a nonzero bulk cosmological constant by replacing 
M 2 with M 2 + £TZ and recovering the Ricci tensor term in Eq. Q . 



1 Generically, s urface terms ar e irrelevant for the Casimir force, but are essential to relate the vacuum energy density and the Casimir 
energy, see |H E3 El Ell] . 

2 We omit the anomaly term since it vanishes in the bulk for odd dimension, and on the brane it can be absorbed in a renormalization 
of the brane tension. 
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2.1. Spectrum 



The Klein-Gordon equation following from the action (0) is separable into radial and dS parts so we introduce the 
mode decomposition <f> = Y] J U p {r) 3^p£m(x,^), where m is a multiple index. The radial equation is 



n+l 



Or 



1 



P 



n 

T 



Up(r) = 



r r' 

while the brane terms can be encoded in the boundary condition 

(d r + ^cs)Kp\ r=ro = 0, 

where symmetry has been imposed and the effective brane mass ^t e ff is given in Eq. J3J. 
The de Sitter part satisfies 



(8) 



(9) 



D n+1 ~(n/2) 2 -p 2 



ypim — o 



Thus one obtains a tower of modes y p e m in dS with masses m^ K = (n/2) 2 + p 2 in units of H. The mass spectrum 
determined by the Schrodinger problem defined by Eqs. JSJ and ©. ft consists of a bound state plus a continuum of 
KK states with p > (jtikk > n l^)- The radial part for the KK modes is of the form 



- r ~n/2 



[A p I ip (Mr) + B p I- ip (Mr)} , 



(10) 



with A p and B p determined by the boundary condition @ and continuum normalization, 2 J drr n 1 Wp K (r)Wp, K (r) = 
S(p-p>). 

The mass of the discrete spectrum is 

m 2 = (n/2) 2 +p 2 <(n/2) 2 , 
and hence pa is pure imaginary. The normalizability implies that its wave function is 

U hs {r) = N d r- n l 2 I_ lpd (Mr) , (11) 
with —ipd > 0. The boundary condition @ determines p d in terms of M and /^ e ff according to 

v I- ipd (Afr ) - Mr I'_ ipd (Mr ) = , (12) 

where we introduced the combination 

v =1- M. (13) 
2 2H y 1 

In the limit Mr a <C 1 and /i e sro <S 1, Eq. (|12|l implies that the mass of the bound state is 



(Hm d ) 2 = nH [i oS + 



M 2 + OU 2 eS , r i eS M 2 ,M i ), 



(14) 



jUeff 




FIG. 2: The shaded area corresponds to the values of /x e ff and M 2 for which the bound state is normalizable (rrid < n/2) and 
non tachyonic, md > 0. The thick (red) line corresponds to the massless case. The plot is for n — 3. 



5 



which agrees with the results of Ref. jToL |2S| . Figure |3 shows the values of M 2 and p c g for which there exists a 
non-tachyonic (~ipd _ n /%) bound state. In this paper, we are mostly interested in the situation when the bound 
state is massless. This happens whenever Eq. (|12|) with —ipd replaced by n/2 holds, that is when v reaches the 
'critical' value 

Mr I' n/2 (Mr ) 

VC = /„ /2 (Mr ) ■ (15) 



2.2. Green's function 



The renormalized Z?-dimensional Green's function can be split into the bound state and KK contributions 

(ren) 



G< ilL = GKK + Gbs . ( 16 ) 



with 

G bs = U bs (r)U bs {r')G^ (dSV (17) 



dp[U™(r)U™(r')] ren GU d3) , (18) 



where G^ dS ^ denotes the Green's function of a field with mass {n/2) 2 +p 2 in n+ 1 dimensional dS space with H = 1. 
It depends on x and x' through the invariant distance in dS, which we call £. Its precise form is given in Appendix 
lAl The 'renormalized' product of the KK mode functions in Eq. (|18|) is 

[U; K (r)U; K (r')] ren ee U; K (r)U; K (r')~U^ Imk (r)U^ mk (r') 

in yK-UMro) - Mr K' . (Mr ) 

Here, U p ' Imk (r) oc Ki p (Mr) is the Minkowski counter part of I|1U|) . This effectively removes the UV divergent contri- 
bution to the Green's function and guarantees the renormalized Green's function (|16|) to be finite in the coincidence 
limit. 

Since Eq. (JTHJ is even in p, Eq. (TBJ can be cast as G KK = dp p™{r)U™(r')] r ™ \ G^ ds) , where 

[W™(r)^ K (r')]^ e " stands for the first term in Eq. (JTJJl only. This can be evaluated summing the residues. From 
Eq. i|A5(l . the poles in G^Xg*. in the upper p plane are at p — i(q + n/2), with q = 0, 1, 2 . . . (see Fig. [3J. From 
Eqs. I|19|) and i|12|) . we see that the KK radial part has a pole at the value of p corresponding to the bound state, 
p = Pd- We shall now show that the residue is related to the bound state wave function as 

2m Res [U™{r)U™{r% en | p<j = -U hs {r)U bs {r 1 ) . (20) 
Using the Wronskian relation K x (z)I' x (z) — K' x (z)I x (z) =l/z and Eq. (|12fl . it is straightforward to show that 
M Res [U™(r)U™(r>)r n \ = -— ( Pd /I- w AMr ) I- iPd (Mr) J ^ (Mr*) 

P — Pd 

The overall constant in the r.h.s. is nothing but the normalization constant in the bound state wave function 119|l . 
up to the sign. Using the Schrodinger equation (JSJ) and integrating by parts, we have 

f r ° dr 

(P 2 -P 2 d) / -I-i Pd (Mr)I- ip (Mr) = I^ p {Mr )Mr l'_ iPi (Mr ) - I- ipd (Mr o )Mr o lL ip (Mr ). 
Jo r 

Setting p = pd after differentiation with respect to p, we find 

1 f r ° dr 1 

-2=2/ - [I- iPd {Mr)f = -I- iPd (Mr ) d p (vI- ip (Mr ) - Mr I'_ ip (Mr )) | , (21) 

+ *d JO r Pd 

where we used Eq. (|12fl . From this and the form of the bound state wave function (flip, it is clear that Eq. I|2U|) holds. 
Equation <|2lTfl implies that no term of the form U bs (r)U bs (r')G p d d S ^ survives in the result. This is 'fortunate' because 
close to r = 0, U hs ~ r ~ n / 2 ~ Wd which is divergent. Thus, Eq. I|2U|I guarantees that (T^) is regular on the light cone. 
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Since only the poles from G^ hg* c ontribute to G^ ren ^ , it can be written as the integral over a contour C that runs 
above pd (see Fig. El- Equation (|A5(I leads to the expression appropriate to see the coincidence limit, 



G(™> = dp [U; K {r)U; K {r% en G« 

J C 



p(dS) 



S (n) _ vK n/2+k+j (Mr ) - Mr Q K' n/2+k+j (Mr ) 



2 „-ir(»±l) ^ n/2+fe+J .(Mr )-Mr /; /2+)b+j (Mro) 

(f_+fc + j) In/2+k+j{Mr)I n/2+k+J (Mr') (-l) k T (n + 2fc + j) fl-cosC 
(rr')"/ 2 j!fc!r(s±i + fc) 



(22) 



and we remind that £ is the invariant distance in dS. Each term comes from the pole at p = i((k + j) + n/2). Setting 
k = 0, we find that the terms with a large j is unsuppressed for r = r' = ro for C = 0. Hence the Green's function in 
the coincidence limit is divergent on the brane. This is the usual UV 'Casimir' divergence near the boundary. Since 
we are interested in the IR behavior, we shall not further comment on this UV divergence. 

The term with k = j = in Eq. (|22|l renders the Green's function globally IR divergent in the limit when the bound 
state is massless, v — > v c (see Eq. (|15|) ") . One can show that this term comes from the homogeneous (£ = 0) mode of 
the bound state. Using Eqs. I|C2() and (|C3|) . the leading behavior of Eq. lll'L'1) in the massless limit nid — > can be 
written as 

G (ren)(X) = ^_ & U b Q * {r)U b Q * {r') + O(m d ) , (23) 



£(ra+l) 



■III 



where U^ s {r) = NqL^ (Mr)/r™/ 2 is the wave function of the bound state for the exactly massless case. 

The divergence (|23f) appears because in the massless limit, the wave function of the homogeneous mode of the 
bound state in the dS invariant vacuum broadens without bou nd Il5j| . It can be removed by considering another 
vacuum with finite width, which implies breaking dS symmetry [14j . Later, we shall take the Allen Follaci vacuum 
[l4|. We will find in Section [3] that in the brane world context, this process removes the global IR divergence, but a 
localized singularity within the bulk remains. 

Let us now examine the behavior of (T^ u ) in the massless limit, wid — * 0. In this limit, the stress tensor is given by 

{T,^ lk ~ g2 2 T M „ [Z^(rKV)] , (24) 

t)(n +1 )m d 

(T tJ ) br ^ ~ - g2 2 [(44 - lKff ~ 2Hi] {U* s | ro ) 2 5{r - r ) % , (25) 

where T^ v is the differential operator given in Eq. 10 and h%j is the induced metric on the brane. For M/fl, Uq s is 
not constant. Therefore Eq. (|24|) explicitly shows the presence of a global IR divergence in (T^ u ) for the dS invariant 
vacuum in the — » limit. For M = 0, Uq s is constant. Hence, the bulk part is finite. However, if £ ^ then 
the surface term (|25|l diverges. 3 Thus, in the limit md — we are forced to consider another quantum state. This is 
the subject of Section We shall mention that there is a possibility to avoid this IR divergence without modifying 
the choice of dS invariant vacuum state. In the present case the divergence is constant and is localized on the brane. 
Hence, it can be absorbed by changing the brane tension. We may therefore have a model in which this singular 
term is appropriately renormalized so as not to diverge in the md — » limit. Of course, such a model is a completely 
different model from the original one without this IR renormalization. 

The massless minimally coupled limit, M = /i e g = £ = 0, is exceptional. Both bulk and brane parts of stress tensor 
(|24|l and 125|) are finite in this limit, though here there is a slight subtlety. The limiting values depend on how we fix 
the ratios among M,/j, e g and £. For example, using Eq. 1)14(1 ■ the surface term is given by 

I- irp \brane ^ U + 2 H 2 (U^ S \ ro ) 2 /I eS + 2gff 

M^o^ - — n (n + 2)H, cS + M 2 5{r ~ ro) ' (26) 

where we used the approximate mass of the bound state, 

m 2 d ~nfiH + 2n(n + 1) £ H 2 H — M 2 , 



3 For M = 0, /i e ff (4> 2 ) is finite because (</> 2 ) ~ l/[i e g, see ('Ea. 1141 . 
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)! i (n/2 +2) 



x i(n/2+l) 



K in/2 




C 



FIG. 3: Integration contour in Eq. 1221 . 



which is valid in the massless minimal coupling limit (see Eq. Q14JI). Hence, in the absence of any fine tuning 
(m^ w max(M 2 ,t;H 2 ,Hn)), it is clear that the contribution (|25(l is not large, even though the Green's function 1231) 
is. Thus, only in the case when the parameters are 'fine tuned' according to Eq. (|f 5(1 the stress tensor l)25|l is large. 
From Eq. ©, if the Green's function is free from 1R divergence, it is clear that the brane stress tensor must be zero 
in the massless minimally coupled case. The direction that reproduces this result is the one along which {Tij) brane 
already vanishes (in the massive case), that is fj, s g = — 2£-ff. Note that this feature is analogous to what happens in 



The bulk stress tensor ()24f) also has a similar but a bit more complicated feature. The operator JJjJ has terms which 
do not manifestly involve a small quantity, such as Af 2 , £ or \i. However, these terms are associated with derivative 
operators. In the limit M 2 — > 0, we can expand the r-dependence of the term with k = j = in Eq. (|22|l as 



Then, the leading term in the above expansion, which is not suppressed by a factor M 2 , vanishes in l|24|) . The 
remaining terms are finite unless the parameters are fine tuned, as in the case of the brane stress tensor. 

Finally, Eq. I|24|l also shows that (T M „) is perfectly regular on the light cone for ma ^ 0. As mentioned before, this 
happens thanks to the KK modes. Note as well that for M > H, Eq. (|24|l is exponentially localized on the brane, 
because the bound state is localized in this case. 



In the preceding section, we have seen that de Sitter (dS) invariant vacuum causes divergence in the limit when the 
bound state is massless. The divergence is caused by the £ = homogeneous mode in the bound state. In this section 
we consider a different choice of vacuum state for this mode aiming at resolving the problem of divergence, following 
the standard methods used in dS space For simplicity, we concentrate on the case of exactly massless bound 

state rrid — that is pd = ni/2, although more general cases would be treated in a similar way. The case with = 
includes not only a massless minimally coupled scalar, M = fj, = £ = 0, but also other fine tuned cases. 



dS space [l5|. 




(27) 



3. 



EXACTLY MASSLESS BOUND STATE 
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3.1. The Green's function 



Here, we should split the Green's function into KK and bound state contributions, G( ren )( 1 ) = G KK ^' + G bs ^ . We 
leave the quantum state for the KK contribution untouched, and change only the contribution from the bound state. 
In the integral representation for Green's function in Q22[l. we have used the integration contour given in Fig. EI This 
choice of contour automatically takes into account the bound state contribution simultaneously. In the present case, 
we consider the contour that runs below the pole at p — pd to exclude the bound state contribution. For md = 0, the 
integrand has a double pole at p = ni/2 because the pole in the radial modes coincides with one of the poles in the 
dS Green's function. Hence the integral with the contour given in Fig. [3] which runs through these merging poles, is 
not well defined. But the integral with the new contour, which picks up the contribution from the KK modes only, is 
well behaved, and it can be cast as 



£jKK(l) 



simple poles ^ d 

hRes (G^)j p ((p - Pd ) [U™{r)U™{r')} [ en ) | 1. (28) 



As before, [U p K (r)U p K (r')~\ r denotes the first term in Eq. (|19|l . In the first term the 'simple poles' mean poles at 
p = i(q + n/2) with q = 1,2,... (see Fig. Namely, it is obtained by removing the term with k = j = from 
Eq. I|22|) . The last two terms are contributions from the double pole at p = ni/2. 

Next, we consider the contribution from bound state. A massless bound state behaves as a massless scalar from the 
viewpoint of n + 1 dimensional dS space. In dS space it is well known that the dS invariant Green's function diverges 
in the massless limit because of the £ = homogeneous mode ■ The usual procedure is to treat separately the 
I = mode from the rest. It is easy to show that (see Appendix lB|l 

G KD = U b %r)U b ^r'){Y / y P e m (x)y;,M')+yAF(x)yiAx')}+ c.c. (29) 

\r)u bs (r'){a p [(p- Pd )G p ds ^] \ Pd ~d p [{p- Pd )yMx)y;M] \ Pd +y A F(x)nAx')}+o. c . , 



u bs 



where y p e m {x) are the positive frequency dS invariant vacuum modes with mass p 2 + (n/2) 2 . The last term in Eq. lj2l3|) 
is the contribution from the homogeneous mode in the appropriate state, which we shall take to be the Allen Follaci 
(AF) vacuum (see Eq. ®). 

At this stage, we note that the second term of Eq. I)28|l cancels the first term of Eq. I|29|l . due to Eq. I|2(JI) . This 
cancelation resembles the one that occurred in the previous case between the KK modes and the bound state. In that 
case, it guaranteed the absence of light cone divergences. In the present case, the terms that cancel in Eqs. <|2~5|) and 
(I29|l are already regular. We show below (see discussion after Eq. \62\i ) that instead the last term in Eq. I|28|l and the 
second term in Eq. I|29|l diverge on the light cone. However, when added up, they render Q ren finite on the light cone 
in odd dimension. (In even dimensions, G ren is finite but its derivatives are not.) 

The fact that the dS invariant Green's function diverges because of the homogeneous mode implies that 
i?es(Gp d5 ^^) pd = liiii p ^ Pd (p — pd)y p ooy p oo +c.c. Using this and Eq. J20J, we can rewrite the total Green's function 
in the more convenient form 

G (re„)(l) =G (l) +Glc + Q AFj {30) 



with 



+ c.c. 



g lc = 2 m d p ([p-pd)y P ^y; m ( P ~ Pd ) [u; K (r)u; K (r% e ^ 

G AF = U bs (r)U bs (r')y AF y'X F +c.c. , (31) 

and G$mvle 1S * ne contribution from the first term in Eq. I|28|l . As is manifest from the expression (|22|) with the 
k = j = term removed, G^ mple is regular in the coincidence limit, everywhere except for the 'Casimir' divergences 
on the brane, and depends on the points x and x' through r, r' and £(ir, x') (the invariant distance in dS between the 
projections of points x and x'), and hence it is dS invariant. The second term, Glc, contains the two contributions 
in Eqs. 128|) and (|29|l that are separately divergent on the light cone mentioned in the previous paragraph. The last 
term, Gaf, encodes the choice of vacuum for the zero mode of the bound state. 
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Let us examine the contribution potentially divergent on the light cone Glc- The derivative of the x dependent 
part is obtained from Eqs. I|B4|) and (|B6|I . The derivative of the radial part follows from Eq. (|19|l and d\I\(z) — 
I\{z) logz + f(z), where f{z) is a regular function. It is easy to see that Glc takes the form 

G LC (x,x') = U bs {r)U bs (r') Re \F(x) + Fix')} + regular, (32) 

n6 ( „ +1 ) 

with 

F(x) = logr + n / / d Xa cosh™ x 2 , (33) 

Jo cosh Xi J -if 

where the indicated regular term depends on r and r' only. The double integral in (|33(l grows linearly with \ for large 
X and eventually blows up on the light cone (The first integral asymptotically grows like e nXl /n, and therefore the 
integrand in the second integral goes to a constant). This is the expected behavior from the massless bound state. 
On the other hand, the KK modes contribute the logr term, which cancels the light cone divergence. To see this, we 
integrate by parts to obtain 

F(x) = log r + log sinh x + f dx 2 ^jj^ , (34) 

Jo cosh Xi J smh x 2 

and now the double integral is bounded. The first two terms are simply log T. Thus, the leading divergence in Glc 
on the light cone cancels between contributions from bound state and from the KK modes, although it still diverges 
logarithmically at infinity. 

This statement has to be qualified for even dimension. In this case, the derivatives of Glc diverge on the light cone 
because of the last term in Eq. I|34l) . To see this, note that the integrand in Eq. 13411 can be expanded in exponentials. 
Then, the integral is a sum of exponential terms except for one, of the form xe _ " x , if n is even. In terms of the null 
coordinates U and V, this is ~ (V/U) n ^ 2 \og(V/U) for x - * +oo (for x ~^ ~°° replace U <-» V). Even though the 
Green's function is regular at V — 0, the stress tensor develops a singularity which behaves like ~ 1/V or 1/U on the 
light cone in four dimensions (n = 2) and like ~ \ogV or \ogU in six dimensions (n = 4) if £ ^ 0. 

For reference, we show the explicit form of F{x) for dimensions 4 and 5, 

2 , 

(35) 

3 , 

where r 2 = T 2 — R 2 . Note that despite appearances, Eq. is regular at R = 0. This is guaranteed since the 
X dependent part of F(x) is related to the dS invariant vacuum modes (see Eq. (|B6ll ). which are regular at R = 
(£ = — iri/2) by construction. The expressions l|35|l are appropriate in the Milne region, and we have already taken 
the real part, which is the relevant part for G^\ 




UlogU + V\og\V\ 

— tor n = 

R 

log (T + r) - (T ; 2 T)T for n = 



The last term Gaf in Eq. (|30|l corresponds to the choice of vacuum for the t = mode. This mode is peculiar 
because it behaves like a free particle rather than an oscillator. The eigenstates of the Hamiltonian, and in particular 
the ground state, are plane waves in field space. However, such states are not normalizable. One can construct well 
defined states as wave packets. The simplest option is a Gaussian packet. This is the Allen Follaci vacuum which 
in fact is a two-parameter family of vacua. Its mode function is given by |14| 



y AF (x) = 



1 ■ r dx ' 

— +i(3-ia / — - 
2a J cosh x 



(36) 



where a > and /3 are the mentioned free (real) parameters. We shall impose that the vacuum is time reversal 
symmetric. This translates into 3^U_f(~ x) = ^af(x)j and implies = 0. Because of the time dependence, it breaks 
dS symmetry. For this vacuum, 



Gaf(x, x') 



ir ' [r]ir ' [lJ] f ^L + 2a 2 Re[F AF (x)FX F (x')} -Im[F AF (x) + F AF (x')} 



(») 



(37) 
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with 



dx cosh n x- 



In order to obtain the analytic continuation of Faf to the Milne region, it is better to write it as an integral over a 
contour with constant T . It is straightforward to see that 



F AF (x) 



dR' 
R 7 ™ 



R'- 



T 



2\ 2 ■ 



It is transparent now that Faf behaves like 1/i?™ -1 near R = 0. It is also clear that it is regular on the light cone 
and is bounded at infinity. Moreover, it gets an imaginary part in the Milne region for odd dimensions (in the Rindlcr 
region it is always real). More specifically, we have 



Faf(x) = 



T 

R ' 

i 

2 



Tt , f.T + T 

— - — In i 

R 2 \ R 



for 
for 



(38) 



3.2. Divergent stress tensor 

Now we discuss the form of the expectation value of the stress tensor for the possible values of M, /i G ff and £ in 
which the bound state is massless. The bulk part of the stress tensor is most conveniently separated into 



v / simple 

(39) 

where (T^ u )q contains the contributions from Glc and G A f , whereas (T^„) simple is the contribution from G^ mple . 
All the IR irregularities are contained in (T^q . From Eqs. (|32|l and Ij37|l . we obtain the bulk part of {T al/ )o as 



jrp \bulk 



1 



T 

- 1 nv 



U bs (r)U bs (r') 



2S'(n) 

-2£ W^ v U bs {r)U bs {r') 



1 

2d 2 
ReF(x 



2a 2 Re[F A F(x)F A _ F (x')]} 
1m Faf (x] 



iS, 



(n+l) 



25" 



(n) 



(40) 



The term proportional to a 2 diverges at R = like 1/R 2n for any value of £. R — corresponds to a timelike axis 
in the bulk passing through the center of symmetry (see Fig. EJ). The last term also diverges like 1/R n+1 for odd 
dimensions, but vanishes for even dimensions. The piece involving Re_F diverges on the light cone for n = 2 or 4. 

Let us begin with the most general case with M =^ 0. In this case U bs (r) is not constant, and the contribution from 
the term inversely proportional to a 2 in Eq. I|40|) does not vanish. This term is analogous to Eq. I|24l) , and it diverges 
globally in the a — > limit. Hence, this state cannot be taken, and one has to content with the AF vacuum with 
some nonzero a. Hence, the term proportional to a 2 is unavoidable. But this is very noticeable since it contains a 
singularity at R — of the form ~ a 2 /R 2n , present even for minimal coupling. The main point is that in the presence 
of a bulk mass, (T^ v ) contains a quite severe bulk singularity even after we get rid of the global IR divergence. There 
is of course the possibility that a different choice of vacua for the KK modes could cancel it out. In this case, it 
seems that the vacuum choice should not be dS invariant, otherwise the singular zero mode contribution that is not 
dS invariant could not be compensated. We leave this issue for future research. 

Now we turn to the case with M — 0. Let us first consider the nonminimal case £ ^ 0. Since M = 0, U hs (r) 
is constant and the bulk stress tensor does not diverge globally in the a — > limit. Divergences in the other a 
independent two terms in (|4UH also vanish in even dimensions with n > 6. However, we cannot avoid the divergence 
in the brane stress tensor in the a — » limit. From Eq. (J7J) and taking into account that ^ ff = (see Eq. d}), we 
have 



brane 



-2^ 2 ) AF ICij S(r - r ) 
£n f (2a)- 2 + a 2 Fl F (x) 2ReF{x) 



»(n) 



'(n+l) 



hij 6(r 



ro) 



(41) 
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where we used that U 2 S — n/2r^. Thus, one is forced to take the AF vacuum again. As a result, the bulk stress tensor 
develops the same singularity at R = as before in the M ^ case. The globally divergent (1/a 2 ) term in Eq. ^IJ is 
proportional to the induced metric. One might wonder whether this effect is physical or not, since it could be simply 
absorbed in the brane tension as before. We think that such a procedure is not justified here because a is a state 
dependent parameter, and renormalization should be done independently of the choice of the quantum state. 

Before we examine the stress tensor in the massless minimal coupling case, let us comment on the relation between 
the AF vacuum and the Garriga Kirsten (GK) vacuum. The latter was introduced in 1 1 ol | for a massless minimally 
coupled scalar in dS. It corresponds to the plane wave state with zero momentum in field space. This is intrinsically 
ill defined, giving a 'constant infinite' contribution to the Green's function. However, the point is that, since it is 
an eigenstate of the Hamiltonian, it does not depend on x so it is dS invariant. Divergence in the Green's function 
can be accepted in the massless minimally coupled case. Our reasoning is as follows. In this case, the action has the 
symmetry <p —> <p+ constant. If we consider it as a 'gauge' symmetry, all the obscrvablcs arc to be constructed from 
derivatives (or differences) of the field. In fact the stress tensor operator JBJ only contains derivatives of the type 
) in this case. Then we will find that the constant contributions in the Green's function are irrelevant. 4 
We need a practical way to compute quantities in this vacuum. If we follow the argument by Kirsten and Garriga, it 
will be given by the limit of AF vacuum 

\0)gk = Urn |0}af • 

From Eq. it follows that the brane term (Tij) brane vanish in the GK vacuum. The first termin Eq. l|4"0)l. which 
is responsible for the 'infinite constant', vanishes for M = 0. While the second term proportional to a 2 also vanishes 
in the a — > limit. We shall note that the last two terms in Eq. (|40|l independent of a are also zero for £ = 0. Hence, 
we are left with (T^) simple, which is manifestly dS invariant, and is finite (aside from the 'Casimir' divergences on 
the brane). Thus, the total stress tensor in the GK vacuum is given by a simple formula presented below in Eq. H5()(l 
with £ = 0. In contrast to the massless minimal coupling limit discussed in the preceding section, here we do not have 
any ambiguity. 



4. ZERO BULK MASS 



In the absence of bulk mass M, the Green's function and the stress tensor can be obtained explicitly. We shall thus 
discuss now this case. 



4.1. Generic mass of the bound state 

First we discuss the case when the bound state is not massless. The other case is postponed until Section l4~2l 
For the bound state, the boundary condition J§J) fixes pd = iv and md is given by Eq. (|14fl with M = 0. Its wave 
function is proportional to r~ lVd ~ n l 2 — r~^ c!f / H , which is constant if = (see Eq. fl3jl ). From the discussion in 
Section l2~Tl the bound state is normalizable and non-tachyonic for n/2 > v > (i.e. < /i c ff < (n/2)H), see Fig. EI 
The normalized radial KK modes are 



OO = ^ rn{l + {v/m (cos bin,) ^(plnr)) , (42) 
and the renormalized product of mode functions in Eq. (|18|) is 

(43) 



27rr n/2 



p + iv , • p — iv . .. ,„ 

(rr') lp + (rr'Y p 

p — iv p + iv 



Proceeding as in Eq. (|22|l . we can explicitly perform the integration over the KK modes and in this case, we obtain 
a simpler expression, 



D« ( XX >)~ 1 y^ j+^+fc + i (-l) fc r(n + 2fc + ,) (rr>\ k+ > ( l-cos( ] 



4 Taking this vacuum is analogous to performing a Gupta-Bleuler quantization 121-1 . 
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where as usual, we denote the (bulk-) massless Green's function by DY^Jx, x'). From Eq. (|44() . it is manifest that 
the Green's function in the dS invariant vacuum state is regular at r = in the coincidence limit. As before, in the 
massless limit v — > n/2, it is singular. However, since the divergent term k = j = is constant in the present case, it 
will not affect the bulk (T M „). 

In order to compute the stress tensor, it is convenient to rewrite Eq. (|44|l in a more compact form. For a conformally 
coupled field, v vanishes, and the Green's function actually takes the simple form 5 



D 



(ren)(l) 



{x,x') 



S (n+i) \r% + {rr' /r ) 2 ~ 2rr' cos ( 



i/2 



(45) 



This expression can be obtained by the method of images. It corresponds to the potential induced by a source of unit 
charge at x' together with an image source located at r\ = r^/r' with a charge q'j = (r'j/ro) n . From the form of 144fl . 
the Green's function for v ^ can be obtained from Eq. (|45(l by applying the integral operator 



£)(ren)(l) 



1 + 4^r n 



ro 



~2v 
r 



D 



(ren)(l) 



(46) 



Borrowing the intuition from the method of images, Eq. (|46|l can be interpreted as the potential induced by the 
image charge mentioned in the previous paragraph together with a string stretching from rg/r' to infinity along the 
radial direction defined by x 1 with a charge line density given by A(r) = 4i/(r/ro) 2l/_1 /r'. To obtain the stress tensor, 
we can first compute the case v = and then apply the same operator as in Eq. I|46|l . The general result with M = 
and rrid ^ is 



(-l) n (n + l) 



„™+2 



*(n+l) r 

(m d¥0){T i ) = (€ _ €e) (-D"(" +1) 



r, 



2n+2 




n 2\n+l 



2F1 [n + 1, - + 1 - v. - + 2 

-2-2I/ V 2 2 



„2n+2/ 2 



(rg + r 2 ) 



4;^ 



(rg - r 2 ) 



2\n+2 



2-2i/ 
,..2 



?2 7* 

2F1 I n + 1, - + 1 - v; - + 2 - v; 



r „ 1 n 
2 F 1 [n + 2,- 



n 

^2 +3 



(47) 



where ^ c = n/4(n + 1). As an aside, we shall note as well that for conformal coupling, (!),„) = even with a 
nonzero boundary mass /i. This is a consequence of conservation and tracelessness of T^ v . 



4.2. Massless bound state 



We consider now the case = 0, that is v = n/2. We can proceed as in Eq. (|30|l and decompose D ren = 
Dsimpie + Glc + G af , where D s %m-pie is the contribution from the simple poles p — i(n/2 + k) with k = 1, 2, . . . in 
Eq. I|28|) . Thus, it is given by the terms in Eq. 144|) with non- vanishing k and j. The integral representation analogous 
to Eq. is now 



D 



(i) 

simple 



2n 



Jr 



ro 



-„ 



D 



(ren) 
u=0 



1 



n5(„+i)f l 



(48) 



where we subtract the constant to remove the j = k = term. The explicit expression in four dimensions was obtained 
in [2l]] 6 , and we shall not reproduce it here. The case of main interest for us is n — 3 and we find, up to an finite 
constant, 



D [ 'L le (x,x') = Re 



simple 



1 



1 



A 3 / 2 



87r 2 r 3 



cos 2C, 
sin 2 C 



2C - {rr'/rl) cos C 
sin 2 C A 1 / 2 



log 



•cosC + A 1 / 2 



, (49) 



5 From Eq. 1441 . the expression for a Dirichlet scalar (v — > 00) is the same with opposite sign. 

6 In Ref. Z2 reflection symmetry is not assumed for the field. 
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where A = 1 + (rr'/r 2 ) 2 — 2(rr' /r 2 /) cos (,, As mentioned above, in the coincidence limit this contribution is regular 
except on the brane. We shall note as well that it grows logarithmically at infinity. 



The contribution (FtivYsilnpie is easily found exploiting the integral representation (J2SJ, as before. The only difference 
between Ij48f) and (|46Jl in the limit v = n/2 (i.e. m,d = 0) is a constant, which does not affect (T M1/ ) in this case. Thus, 
{T^) b ^impie reduces to (1471 with v = n/2. This gives a simple expression in terms of elementary functions. In four 
dimensions it is given in |2l| (for £ = 0). In five dimensions, we obtain 

9 £"£c f rl , r 4 -3r 2 r 2 +3r 4 



( r pr \ bulk | _ 

\ r/simple ~ 327r 2 r 3 ^ ^2 _ r 2)4 ^ ( f 2 _ r 2)3 

lT i,bulk = 9 ^c| 6 ^ + r| r 2 r 4 - 4r 2 r 2 + 6r 4 r 4 - 3r 2 r 2 + 3r 4 ) 

\ i/«mpie 327r 2 r 3 | r ( f 2 _ r 2)5 + 2 (r 2 - r 2 ) 4 (r 2 - r 2 ) 3 J J ' ' 1 j 

which is dS invariant and regular everywhere except on the brane. 



5. PERTURBATIONS 



We shall now discuss the form of the field perturbations, focusing on the bulk massless case with n = 3, since we 
have obtained a closed form expression for the Green's function and the stress tensor in the preceding section. We 
begin by the case with generic coupling, and we consider (4> 2 (x)) in the Allen Follaci vacuum. The massless minimally 
coupled case in the de Sitter (dS) invariant Garriga Kirsten vacuum requires a different treatment, which will be 
discussed in Section l5~2l 



5.1. Generic coupling 

The renormalized expectation value of <fi 2 {x) is given by 

with Grin given by l|3Q(l . From Eqs. (|37(l and (|38|l . it diverges at R = because of Gaf- It is also clear that 
this contribution is bounded at (null) infinity. From Eq. I|32|l . the Glc contribution is regular on the light cone. 
Equation 149|) shows that the KK contribution D s i mp i e diverges on the brane. In the bulk, both D S i mp i e and Glc 
grow logarithmically at infinity. However, the growing terms cancel out. Indeed, Eq. (|49|l shows that as long as x is 
not on the brane, D simp i e at coincident points is dominated by the logarithmic term. Therefore we have 



1 



^l vle (x,x) ~ -—^3 log 1 - 3 , (51) 



2 - simpler ■»> 16 ^2 r 3 



r 



,2 



'0 

in the limit x — > oo. As for Glc, since the wave function of the bound state in the M = case is W 2 S = n/2r^, we 
find 

\G LC {x, x) = Re[ log (T -ir)]+ 0(1) , (52) 

where we used Eq. H35(l . It is clear that the logarithmic terms cancel and as a result {4> 2 {x)Yaf lii bounded at infinity. 
This is expected, because the bulk is flat. Intuitively, in four dimensional dS space, the perturbations grow because 
when the modes are stretched to a super-horizon scale, they freeze out. Since modes of ever smaller scales continuously 
being stretched, they pile up at a constant rate Since this effect is due to the local curvature of the spacetime, 
it should not happen in a flat bulk. 

Accordingly, on the brane we recover the same behavior as in de Sitter space. Indeed, restricting (|52|l on the brane, 

we obtain G£q ~ x- We have mentioned before that D S i mp i e (x, x) is UV divergent on the brane. Since this happens 
because point splitting regularization used here does not operate on the brane, this object needs UV regularization 
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and renormalization. This can be done in a variety of ways, e.g. with dimensional regularization 7 , introducing a 
finite brane thickness, smearing the field etc. The point is that the renormalized value must be a constant simply 
because D simp i e (x,x) is dS invariant and is a function of x only. Thus, Eqs. (|5*T|l and (|5*2f imply that on the brane, 
{(j) 2 ) grows in time, as in dS space. 



5.2. Massless minimal coupling 

Now, we shall see essentially the same features arise for a massless minimally coupled scalar in the GK vacuum, 
in which case everything can be obtained in a dS invariant way |15| . Because of the shift symmetry (j) — > 0+const. , 
((j) 2 ) is not an observable in this case. Still, it is possible to define a 'shift-invariant' notion for the field perturbations. 
Following [15j , one introduces the correlator 

g{x,y) = ([<t>{x)-cf>{y)] 2 ) GK , (53) 

which can be thought of as the combination [G^(x,x) + G^(y : y) — 2G^ (x, yj\ /2. Since the first two terms are 
UV divergent, we shall rather consider 



G ren (x,y) = l 



G ren(1) (a;, x) + G ren{1) (y, y) - 2G ren(1) (x, y)] . (54) 



The main point is that all the terms of the form f(x) + f(y) in G' re "( 1 ) cancel out in the combination Q ren (x,y). All 
the terms in Eqs. (|32[) and H37[) that do not vanish in the limit a — » are of this form because U hs is constant in 
the massless minimally coupled case. Thus, this correlator is well defined for the GK vacuum [15j and we can readily 
write 

g ren (x,y) = \ [D { s ]l ple (x,x) + D^ mple (y,y) ~ 2DlH ple (x,y)] , (55) 

with D s impie{x , y) given by (|49|l . 8 Thus, the behavior of the perturbations in the GK vacuum for x and y distant is 
thus described by the asymptotic behavior of D S i mp i e (x,x) and D S i mp i e (x,y). The former is summarized in Eq. I|51|l . 
As mentioned before, for finite x, D S i mp i e (x,x) is regular everywhere except on the brane (in which case yi = y). 

Before describing the asymptotic form of D S i mp i e (x,y), we need to discuss the singularities that it contains. The 
combination that we called A in Eq. i|49|) is 

A(x,y) = |x- y/ | 2 /|y/| 2 , 

where \x\ 2 = 'q^x^x 1 ' and yf = (fo/|2/| 2 )y M i s the image of x (see comments around Eq. (|46f) and Fig.JSJl. Note that 
when y is in one of the Milne regions, then yi is in the other one. Rather than an 'image charge', it represents the 
point where the light cone focuses, (see Fig. 01(b)). For x ^ y and both finite, D^ mple (x,y) has singularities in two 
types of situations. One is when x is on the light cone of yz (then, A = 0). The other case arise when the argument 
of the logarithmic term in Eq. (I49() vanishes. It is convenient to rewrite this term as 



8ttM 



log 



(\x - yi\ + \yi\) 2 - \x 



(56) 



The argument vanishes when they are aligned with respect to the origin and \x\ 2 > \yi\ 2 . This condition defines an 
'image string' stretching from yi to infinity. Hence, this singularity occurs only when two points are in different Milne 
regions. This is consistent with the interpretation of Eq. (|46|l . that the Green's function can be constructed with a 



7 It is illustrative to consider the massless conformal case. The Green's function 1451 for x and x 1 on the brane is simply D JJ!q (x,x') ~ 
l/\x — x'\ n / 2 . If n is negative, this is clearly zero in the coincidence limit, and hence the continuation to n = 3 is zero as well. In the 
non-conformal case, we can use the integral form H46I to compute (0 2 )5(r — ro) in arbitrary dimension. Upon continuation to n = 3, 
one obtains a pole ~ l/(n — 3) and a nonzero finite part. The pole can be absorbed in the brane tension and {<f> 2 (ro)) ren is a finite 
constant. 

8 Note that 1551 is completely regular on the light cone. 
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FIG. 4: The singularities in D simp i e (x,y) are located on the light cone emanating from yi, the image point of y. Also when x 
is on the light wedge emanating from 'image string' stretching from yi to infinity. Here, we plot the image point (smaller filled 
circle) with the image string attached to it. In (a), y is spacelike, in (b) null, in (c) timelike and in (d) at null infinity. 



mirror image and a linear charge distribution over such a string, see Fig. 0] 9 These two situations correspond to the 
coincidence of x with the image of y or its 'light cone'. Therefore all of them are of UV nature. 1 

The asymptotic form of D s imple{% , y) for distant points follows from Eq. (|56|l because then, the inverse powers of 
A in (|49l) are finite. Taking y fixed and x — > oo (with x not on the light cone from yi nor on the image string nor on 
its light wedge), one obtains 

DiLpiefa v) - —^s log H ■ (57) 

However, since D S i mp i e (x,x) grows twice as fast (see Eq. (I51ll ). the combination (|54() is bounded in the bulk. 

We can consider as well both x and y approaching null infinity in the bulk. In this case, the image yi approaches 
the light cone, as illustrated in Fig. 0] (d). Then the logarithm in Eq. f57)| is ~ log (\x — yi\/\yi\ 2 ) — log (\x\ \y\), 
hence the combination Q ren is bounded again. 

On the brane, the situation is very different. As mentioned before, D S i mp i e (x, x) is divergent but when properly 
renormalized, it is simply a constant because of dS symmetry. Then, Q ren (x,y) behaves like D S i mp i e (x,y) and from 
Eq. I|56|l for large separation and \x\ = \y\ = rg, one obtains 

^«(z,y)~ ^ \0g\x-y\ . 

Since \x — y\ 2 = 2rg(l — cos £), this corresponds to the linear growth with the invariant dS time interval <^(x, y) between 
x and y [l5|. 

Finally, we shall mention that the generalization of Eq. (|53[) , 

4>{x) cj)(y) 



U bs {x) U bs {y) 



In four dimensions the logarithmic term is log \\x — J//| 2 /|3//| 2 | ■ Thus, this string singularity does not appear, despite the same 

interpretation of Eq. 1461 holds. 

As an aside, let us comment on the singularities present in the bound state contribution. Specifically, in the term dp \(p — p ( j)G^ )— 1 J 

prese nt in Eq . 1291 . We remind that this term is canceled by the KK contribution 1281 . Hence, it plays no role in the Green's function. 
From Il4lll r | we know that for x and y timelike related and far apart, then f is large and pure imaginary, and it behaves like 

8p\<P- PilK- 1 } ~ log 1 1 - cos C| = log I (\x\-\xf-\x-xi\ 2 I . 

This presents divergences on the light cone and whenever x and yj (or equivalently, y) are aligned. This condition reduces to the 
coincidence limit on the brane, but in the bulk one does not expect singularities to appear at all these points. The above mentioned 
cancelation guarantees that these unphysical singularities are not present in G ren (x,y) once we include the KK contribution. 
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is completely regular in the GK vacuum when and its asymptotic behaviour at infinity parallels that of l|53|l . 



AF 




FIG. 5: Divergences in (<ft 2 ) ren = G(x,x) ren and {T (J „) ren when the bound state is massless. For the AF vacuum, both the 
Green's function and the stress tensor diverge at R = 0, indicated with a thick dashed line. In addition, if £ 7^ the stress 
tensor diverges on the light cone like 1/U or logarithmically in four and six dimensions respectively (thin dashed line). Besides, 
we also have the ultraviolet 'Casimir' divergence on the brane, represented by the plain thick line. In the Garriga Kirsten 
vacuum (for the massless minimally coupled scalar), {T M „) ren presents only this Casimir- type divergence on the brane. 



6. CONCLUSIONS 

Our main results can be summarized as follows. In analogy with what happens in de Sitter (dS) space [14J. Il5j. 
scalar fields with a massless bound state in the spectrum do not have a well defined dS invariant vacuum, except for 
the massless minimally coupled case. (The case of vanishing bulk mass with non-vanishing curvature coupling has a 
little subtlety, though.) The Green's function and the v.e.v. of the energy momentum tensor diverges everywhere. 
The simplest alternative from the analogy to dS case is to take the Allen Follaci vacuum. However, in this vacuum, 
divergences in the stress tensor are not removed completely within the bulk. Figure [S] illustrates the location of the 
IR singularities in the Green's function and the stress tensor for the AF vacuum. It remains to clarify whether it 
is possible or not to avoid these singularities by choosing a vacuum for the KK modes other than the dS invariant 
vacuum. 

When the bound state is very light (but not exactly massless) because M, fj, and/or £ are fine tuned according to 
Eq. (|15fl . then the stress tensor in dS invariant vacuum takes the form of (|24l) . The stress tensor in this case is smooth, 
but it becomes very large. Hence, even when the bound state mass is not exactly zero, the dS invariant vacuum looks 
problematic because of large back reaction. Note that the situation here is different from the usual dS case in two 
respects. In dS space the large v.e.v. in the stress tensor for the dS invariant vacuum is a constant proportional to 
the metric. Hence, it might be absorbed by IR renormalization of the cosmological constant. In our case, the stress 
tensor given by Eq. I|24l) is a nonlocal expression and cannot be 'renormalized away'. On the other hand, if one does 
not want to make any IR renormalization, in the dS case one can take the AF vacuum, and {T^} stays regular. In the 
brane world, we do not know the prescription how to remove this large v.e.v. by changing the vacuum state. Choosing 
non dS invariant vacuum will lead to not only the mentioned divergence in the bulk, but also a new singularity on the 
light cone, when the bound state mass mj is not exactly zero. In this case, the radial function for the bound state 
behaves like oc r~ rnd/n near the light cone at r = 0. Hence if you single out the contribution to the Green's function 
from the bound state, it is singular and its derivatives diverges at r = 0. Hence, as far as we restrict the change of 
quantum state to the bound state, we will not be able to remove the large v.e.v. of stress tensor without spoiling its 
regularity. 

A light bound state is compatible with a well behaved and not large stress tensor only in a situation 'close' to the 
massless minimally coupled case. More precisely, here we consider the cases that all of the bulk mass M, the brane 
mass \x and the nonminimal coupling £ are small (see Eq. ()14|l ). This corresponds to having a light bound state 
without accidental cancellations. Namely, the squared bound state mass m 2 , is of order of the largest among M 2 , H/j, 
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and -ff 2 £, where H is the Hubble constant on the brane. In this case, the dangerous terms proportional to m^ 2 are 
always associated with some small factor M 2 , /i or £, and therefore none of them becomes large. 

The case with M 2 ks H/i <C H 2 ^ has a little subtlety. In this case, the large v.e.v. in (T^ v ) appears only in the 
brane part, and it is constant proportional to the induce metric. Hence, we might be able to consider a model with 
an appropriate IR renormalization. Only in such a modified model, the stress tensor in the dS invariant vacuum can 
escape from appearance of a large v.e.v.. 

Application to the bulk inflaton type models 0, El is a part of motivation of the present study. In these models 
there must be a light bound state of a bulk scalar field. In order to explain the smallness of the bound state mass 
it will be natural to assume that it is due to smallness of all the bulk and brane parameters without fine tuning. 
Therefore we will not have to seriously worry about the backreaction of the inflaton in the context of bulk inflaton 
type models. 

We have a few words to add on the massless minimally coupled case. If we consider this case as a limiting situation 
close to the massless minimally coupled case, the results depend on how we fix the ratio amongst M, H^i and H 2 £, 
and hence there remains ambiguity. However, this limiting case has the shift symmetry — > 0+constant. If this shift 
symmetry is one of the symmetries that are to be gauged, there is no ambiguity because the problematic homogeneous 
mode does not exist in the theory from the beginning. In this setup, undifferentiated 4> is not an observable. In fact, 
(Tfiv) automatically does not contain undifferentiated <f>. Hence, (T^ v ) is unambiguously defined although the Green's 
function is not well defined. We can compute (T^ v ) in this model by applying the idea of the Garriga Kirsten vacuum 
(equivalent to a limiting case of AF vacuum), and we confirmed it dS invariant and regular as is expected. 

We have also discussed the form of the field perturbations (0 2 ) when the bound state is massless. The main point 
is that on the brane (4> 2 ) grows linearly with dS time \ while in the bulk it is bounded, as expected since the bulk is 
flat. Aside from this, we derived in closed form the (T^) for a generic field with zero bulk mass, see Eq. 1)47(1. 

The same discussion applies in the RSII model with little modifications, which mainly comes from the fact that the 
Ricci tensor term in the bulk stress tensor ~ ^((f) 2 )^^ is not zero in the RSII model. Thus, the bulk part of (T^) is 
finite in the limit of massless bound state only for massless minimal coupling. 
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APPENDIX A: MASSIVE GREEN'S FUNCTION IN DE SITTER 



In this Appendix we obtain the form of the Green's function for a massive field propagating in D — 1 = n + 1 
dimensional de Sitter (dS) space. The metrics for dS space and its Euclidean version are 



dS 2 n+1) = -dx 2 + cosh xdftfn) = d X 2 E + sin 2 x E dfl 2 n) , 



(Al) 



where dfl 2 n ^ is the metric of a unit n dimensional sphere. The Euclidean time is given by x. E — *X + t/2. The 
Euclidean version of the Hadamard function in dS space can be found from the equation 



d 2 +ncotx E d x „ - [p 2 + — - 



G<g s \x,x') = -6^ +1 \x-x'). 



(A2) 



From the symmetry, we can choose x' to be at the pole so that Gp depends on x E only. In terms of F = 
(sin(x E -7^))( ^l - 1 )/ 2 G^ S) (a;,x , ), this becomes 

(1 - w 2 )dl - 2wd w - 

where w = cos(x E — it), and this is solved to give 



(n- l) 2 
4(1 -w 2 ) 



F = 0, 



(A3) 



(n-l)wi _ I 71 + 1 

F = Ne-^^T 



P I (cos(x E - tt)). 

ty 2 



(A4) 
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Hence, the Green's function will be given by the replacement of x E by the proper distance between the points x and 

x' in dS space, which we call £(ir, x'). G p dS ^ is guaranteed to be regular at £ — > tt. To see the behaviour in the £ — > 
limit, the alternative expression 

is relevant, and here 



N = , M 2 2 j , N . (A6) 



r s-ip r s + ip 



In the coincidence limit, the Green's function must behave like g9^ + ^ w l/(n — l)S( n )C n 1 , where is the 

area of an n-dimensional unit sphere. The limiting behaviour is controlled by the first term. Hence we have ./V = 
2-(»-D/V(n-l)5 (n) . 



APPENDIX B: MASSLESS GREEN'S FUNCTION IN DS 

Here, we compute the Green's function for the massless scalar in dS. In the massless limit (p — > in/2), the de 
Sitter (dS) invariant Green's function diverges because of the contribution from the £ = mode, 3-pOO- The idea is to 
construct a modified Green's function by substituting this mode by another one that is finite in the limit p — > in/2. 
In other words, 

r (dS)(+) _ r {+) ™,» 
Lr (m=0) ~ °(£>0) + yy ' 

where y corresponds to the £ = mode and 

£>0 



where 3Vm are the positive frequency dS invariant vacuum modes. The latter can be obtained as follows. Since 
G p dS ^ diverges because of the £ = mode, the divergent term in the Laurent expansions of both Gp dS ^ and J^ooiVpoo 
coincide. We show below that this series contains a simple pole in p — in/2. Hence, we can write 



G (dS)(+) _ y y 



pOO^pOO 



On 



(Bl) 



The explicit expression for the first term in the r.h.s. is unnecessary because it is canceled by the KK contribution 
}. The dS invariant vacuum mode for £ = with — ip < n/2 is 



3>p00 = C p ^ 2 ^ P> , (i sinhx) 

(cosh X )~ 5 



(B2) 



and the normalization constant has been chosen so that limp^j„/ 2 y p oo/G p = 1. The Klein Gordon normalization 
requires 



r(fr + f)r(-» P +g) 
2 "r(n±i) 2 5 (n) 



(B3) 



Expanding around p = in/2 we find 



|C | 2 = 1 1 

inS {n+ i) p - in/2 



+ 0[(p-in/2)°] 



(B4) 
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where S( n +i) is the area of an n + 1 dimensional sphere of unit radius. The behaviour of y p oo near p = in/2, is most 
easily found using the equation that it solves, 



d 2 + n tanh 



P 



y 



pOO 







(B5) 



Bearing in mind that the positive frequency function for the dS invariant vacuum is determined by the regularity 
when the function is continued to the Euclidean region on the side that contains x — —iri/2, one finds 



dXi 



1 



cosh™ Xi J- 



Xi 



i/2 



d X2 cosh" X2 +o((p-in/2) 2 ) 



(B6) 



From this result, the last term in Eq. I)B1|) can be readily evaluated, and Eq. (j32(l follows. 



APPENDIX C: DIVERGENCE IN THE GREEN'S FUNCTION AND THE I = MODE 



Here, we show that the IR divergence in the Green's function ill'L'l) or l|44(l is due to the homogeneous £ = mode 
of the bound state. Using Eqs. JUJ, l|21|l . 1B4| and l|B6|l . the contribution from the I = mode of the bound state 
for pd close to in/ 2 is found to be 

G ( £ 0<Pd = U bs (r)U bs (r') ^oo(x)^oo(x') + cc. a 2U bs (r)U bs (r') \C Pd \ 2 [1 + ( Pd - in/2)] 

1 1 I n/2 (Mr)I n/2 (Mr>) / (rrT /2 _ 

h . . . (^1) 



S( n+ i) Pd-in/2 I n/2 [Mr ) d p [y c I- iv {Mr ) - Mr I'_ ip (Mr )) 



p—in/2 



where the dots denote higher order in pd — in/2. In the same limit, v is close to v c . The j = k = term in (|22|) is 



(ren)(1) = nT (n) uK n/2 (Mr ) - Mr K' n/2 (Mr ) I n/2 (Mr)I n/2 (Mr') 

3=Q - k=Q 2«r(2±i) 2 5 (n) vI n/2 (Mr ) - Mr r /2 (Mr ) {rr'Y /2 

1 v c K n/2 {Mr ) - Mr K' n/2 (Mr ) + (u - u c ) K n/2 (Mr ) I n/2 (Mr)I n/2 (Mr') 



S {n+ i) I n/2 {Mr ){v-v c ) (rr>r /2 

+ O [{v - v c f] (C2) 



1 I n/2 (Mr)I n/2 (Mr')/(rr') n / 2 



S(n+1) J „/ 2 ( Mr 0) ( V - V c) 

where in the second line we used Eq. (|15fl and the Wronskian relation K n / 2 ( z )I' n / 2 ( z ) — ^-' n / 2 { z )^ n /2^ z ) — Using 
Eq. i(T2|) . it is easy to see that 

n I / .n\ d p [p c I- ip -Mr lL ip (Mr )]\ /2( n , 
v-»c = dvA Pd = m /2 (P* * 2 J + ' ' ■ = I^(M^) ( Pd ~ l V + • • • ' (C3) 

so <|C1|1 and (|C2|I agree in the limit pd in/2. Note that since in this limit this is the dominant contribution and 
— in(pd — in/2), the total Green's function l|22|) can be rewritten in the simple form 

Ggio( ^ )K » !£^) +0( O, ( «, 

J(n+1) ^0 m d 

where U^ s {r) = N Q I n/2 {r) /r n / 2 is the wave function of the bound state for the exactly massless case. 
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